Ground State Properties of Anderson Impurity in a Gapless Host 
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00 ■ Using the Bethe ansatz method, we study the ground state properties of a U — > oo Anderson 

\ impurity in a "gapless" host, where a density of band states vanishes at the Fermi level ep as 

|e — €f\- As in metals, the impurity spin is proven to be screened at arbitrary parameters of the 
system. However, the impurity occupancy as a function of the bare impurity energy is shown to 
acquire novel qualitative features which demonstrate a nonuniversal behavior of the system. The 
latter explains why the Kondo screening is absent (or exists only at quite a large electron-impurity 
coupling) in earlier studies based on scaling arguments. 
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The physics of "gapless" dilute magnetic alloys, where an effective density of band electron states varies near the 
Fermi level as |e — €f\ t , r > 0, has been attracting a significant theoretical interest. Using poor-man's scaling for the 
spin-i Kondo model, Withoff and Fradkin Q| have predicted that the Kondo effect in gapless systems takes place 
only if an effective electron-impurity coupling exceeds some critical value; otherwise, the impurity decouples from 
the electron band. Numerical renormalization group (RG) calculations, large-A^ studies, and quantum Monte Carlo 
simulations have confirmed this prediction and revealed a number of additional intriguing features of the physics 
I , of magnetic impurities in unconventional Fermi systems. 

In a conventional metallic system with (i) a linear dispersion of electrons near the Fermi level, e(k) = vp(k — kp), 
and (ii) an energy independent hybridization, basic "impurity" models are exactly solved by the Bethe ansatz (BA) 
|p[--|ll| . It has recently been shown (l2) also that integrability of the U ~ * oo nondegenerate and degenerate Anderson 
models is not destroyed by a nonlinear dispersion of particles and an energy dependent hybridization, but it becomes 
only hidden |]l3| . The approach developed has allowed us to study |Q the thermodynamic and ground state properties 
of an Anderson impurity embedded in a BCS superconductor, and can be used to obtain an exact solution of the 
■^j- ■ Kondo problem in other unconventional systems. 

In this Letter, we report an exact BA solution of a model describing a U — > oo Anderson impurity embedded in a 
gapless host. The model is diagonalized by BA at the arbitrary density of band states p(e) and hybridization t(e). In 
the RG approach, the physics of the system is assumed to be governed only by an effective electron-impurity coupling 
r(e) = p(e)t 2 (e) rather than by separate forms of p(e) and t(e). However, to derive thermodynamic BA equations, one 
has to specify separate forms of an effective electron-impurity coupling and an inverse dispersion of band states k(e). 
""j""^ ■ While an effective coupling determines the electron-impurity and effective electron-electron scattering amplitudes, an 
inverse dispersion accounts for the spatial behavior of electron wave functions, and naturally enters BA equations via 
{~j periodic boundary conditions imposed on eigenfunctions of the system. The physics of the system is thus governed 
by both an effective electron-impurity coupling r(e) and an inverse dispersion of band electrons k(e) p"5| . 

Here, we treat the case of an energy independent hybridization, t(e) = t = constant, so that an energy dependence 
of an effective coupling T(e) = 2Tp(e), where 2T = t 2 , is determined only by a nonlinear band dispersion. We assume 
O ■ also a simple form for the density of states of a gapless host, 
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where the energy e is taken relative to the Fermi value tp. The parameter (3 characterizes the size of region with an 
unconventional behavior of p(e). At /3 = 0, the model reduces to the metallic Anderson model. If (3 exceeds essentially 
a band half width D, (3 ^> D, one obtains the density of states p(e) ~ |e| r . To derive thermodynamic BA equations 
one has to fix the power r in Eq. (1). The magnitude of r is a key factor in determining the spectrum of the system 
in terms of Bethe excitations. In this Letter, we focus on the simplest case r = 1, which is however of particular 
physical interest [|l]-[jj. 

In the BA technique, the spectrum of elementary excitations of the metallic U — > oo Anderson model is described in 
terms of unpaired charge excitations, charge complexes, and spin excitations including bound spin complexes p[-pT|. 
Since the ground state of the system is composed only of charge complexes (a charge complex contains one spin 
wave and two charge excitations) carrying no spin, the Kondo effect takes place: the impurity spin vanishes at zero 
temperature. 
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In the gapless model described above, the structure of the spectrum is shown to preserve basic characteristic features 
of the metallic version. In particular, the ground state of the system is still composed only of charge complexes carrying 
no spin. Therefore, as in metals, the Kondo screening of the impurity spin takes place at arbitrary parameters of the 
system, that contradicts dramatically to results of earlier studies 0-0]. 

However, the behavior of the impurity occupancy n d as a function of the bare impurity level energy ea is drastically 
changed compared to a metal host. At positive values of ed, the impurity occupancy is still given by the standard 
formulae |^|, where the renormalized impurity level energy e* d = e c i + ^ In — + P contains now the parameter (3. The 
mixed-valence regime shrinks: the impurity occupancy quickly grows from rid ~ at e d = to rid — 1 (precisely!) at 
Ed < — r 2 /4/3. In the Kondo (local-moment) regime, rid = 1 and it does not depend on e d - Only in the empty-level 
regime (e<j > 0) the impurity occupancy is a universal function of the renormalized impurity energy t* d rather than a 
function of the bare parameters of the model. 

The behavior of the impurity occupancy in the mixed-valence and Kondo regimes is not universal that manifests 
nonuniversal properties of a gapless system in contrast to a metallic one. This explains why the Kondo screening is 
absent (or exists only at quite a large electron- impurity coupling) in earlier studies based on scaling arguments. 

We start with the Hamiltonian of the nondegenerate Anderson model rewritten in terms of the Fermi operators 
c\{e) (c CT (e)) which create (annihilate) an electron with spin a —1, J, in an s-wave state of energy e, 



H = H c + H d + H h . (2a) 



Here 



H - = T,I ^ £ 4(e)c CT (e) (2b) 
H d = e d dUa + Ud\d^d\di (2c) 

(7 

H h =Y,I D deVmicUW* + 4c a (e)}. (2d) 

Z J-D 



are the conduction band, impurity and hybridization terms, respectively. All notation in Eqs. (2) are standard. An 
electron localized in an impurity orbital with the energy e d is described by the Fermi operators d a . The electron 
energies and momenta are taken relative to the Fermi values, which are set to be equal to zero. The integration over 
the energy variable e is restricted by the band half width D. In what follows, we assume that D is the largest parameter 
on the energy scale, D — > oo. In the energy representation, the effective particle-impurity coupling T(e) = p(e)t 2 (e) 
combines the density of band states, p(e) = dk/de(k), and the energy dependent hybridization i(e). 

In the limit of a large Coulomb repulsion in an impurity orbital, U 3> -D, eigenvalues of the model (2) with the 
arbitrary inverse dispersion k(e) and effective coupling T(e) are found from the following BA equations [ jl2j : 

t fc,z A' ~e d /2r- i/2 = -p4 hj -\ a - i/2 

hj - e d /2T + i/2 J-l hj - X a + i/2 [ ' 

, K-hj+i/2 11 X a -\p+i 1 ' 

j = l ■> f3=l M 

where N is the total number of electrons on an interval of size L and M is the number of electrons with spin "down" . 
The eigenenergy E and the z component of total spin of the system are found to be 

E = 2Tj2"j, S z = --M. (3c) 

3=1 

In Eqs. (3), lu — e/2T is a dimensionless energy, and kj = k(uij) are charge excitation momenta. From Eq. (1), one 
easily obtains 

k_ _ J w + 51n(l- f), oj < 



2r \ w-ain(l + $), cj>0, (4a) 

where S = (3/2T. The BA equations (3) are quite similar to those in the conventional Anderson model p|-|lT|. The 
only but very essential difference is a nonlinear energy dependence of momenta and charge "rapidities" hj — h(u)j). 
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At arbitrary p(e) and i(e), the rapidity h(e) 
p(e) is denned as in Eq. (1) with r = 1, 



(e — ed)/r(e). In our model of a gapless host, where t 2 (e) = 2T and 



2r' 



As in the metallic Anderson model, in the thermodynamic limit spin rapidities A c 
bound spin complexes of size n, 



(4b) 



1, . . . , M are grouped into 



AL nJ) =A + J(n+l-2j), J 



1. 



(5) 



Apart from unpaired charge excitations with real rapidities hj, the system spectrum contains also charge complexes, 
in which two charge excitations with complex rapidities bound to a spin wave with a real rapidity A Q , 



h^(to) = X a ± 



(6) 



In the absence of an impurity term [the second term in 1. h. s. of Eq. (3a)], BA equations of "impurity" models 
describe a free host in terms of interacting Bethe particles with an arbitrary rapidity h(u>). Introducing an impurity 
fixes an expression for h(uj), and fixes thus the spectrum of Bethe excitations of a host. For instance, in the theory 
of metallic magnetic alloys, the finite-?/ and U — > oo Anderson impurities require different descriptions of the same 
host with the different spectra of Bethe excitations. In a gapless system, the impurity energy e<j is involved not 
only in the impurity term but also in the expression for the rapidity h(u>). The magnitude of td thus determines a 
particle- impurity scattering phase and dictates also an appropriate choice of the Bethe spectrum of a host. 

It is instructive to start our analysis of the ground state properties of the system with the simplest case e<j = 0. 
The complex energies of charge excitations of a charge complex are then found to be 



, X + 5± 5, A< -6 
"±(A)H X -S±l \>6 



(7) 



In contrast to the metallic case, the spectrum of charge complexes contains thus the gap of size 25. The existence 
of the gap should lead to essential changes in the thermodynamics of the system. However, it can be shown from 
the thermodynamic Bethe ansatz equations that the renormalized energies of unpaired charge excitations and spin 
complexes are still positive in the absence of an external magnetic field. Therefore, as in the conventional Anderson 
model, the ground state of the system is composed of charge complexes with negative renormalized energies only. The 
upper edge of filled states, Qg, is now given by Qg = Qa — 5, where Qa — — sr-lnS corresponds to the metallic 
Anderson model. The occupancy of the impurity level is also given by the standard B A formulae M , where in the 



expression for the renormalized impurity level energy e* d one needs only to replace Qa by Qg to get e* d 
Thus, the impurity occupancy at = is essentially decreased compared to the metallic case. 
At e c i 0, the complex energies of charge excitations of a complex are found from the equation 
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(8) 



Since we study here only the ground state of the system, we may restrict our consideration to the solutions with the 
negative real part of energies, Rew < 0. Then, a solution cj±(X) — x(X) ± iy(X) is given by 



where 



i 

sgn/i 



M 2 - b+y/(n*-b)* + n 2 



v(p) 



1/2 



V2 



i 2 + b+ ^{^-by + , 



1/2 



(9a) 

(9b) 
(9c) 
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and /J, = X + 6. The behavior of this solution is governed by the parameter b 

0, the function a;(A) is 

negative only at A < —5. The ground state of the system is still composed of charge complexes filling all the states 
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from A = —D/2T to A = Qg- The impurity occupancy rid is governed by the well known formulae with insignificant 
corrections related to small deviations of the function a;(A) from the linear behavior. 

Significant changes in the behavior of the system occur at negative ed- At negative b (ed/2T < —1/165), the function 
x(X) is negative at all A £ (—00, 00). In contrast both to the metallic model and to the gapless system with positive 
ed, the bare energy of charge complexes, £o(A) = 4Fx(X), is now monotonically increasing negative function at all A. 
Therefore, in the ground state of the system charge complexes fill out all allowed states on the A axis. The density of 
states of charge complexes er(A) is found from the continuous limit of Eqs. (3), 

' ' /,/!A! \a(\-£) = v{\) = [°° dX'a(X-X')a(X'), (10) 



2tt dX L 2T' 

where q(X) = fc-(A) + fc+(A) is the momentum of the charge complexes, and a(x) = [ir(x 2 + l)] -1 . As usual, the 
function er(A) is divided into the host and impurity parts, c(A) = ct/j(A) + L~ 1 <Ji(X). The occupancy of the impurity 
level, rid, is then given by 



n d = 2 / d\ai{X), (11) 
J ~oc 

where the impurity density of states is found from the equation 

poo 

a(X - ^) = ai(X) + J dX'a(X — A')fj(A'). (12) 

Since unpaired charge excitations and spin complexes are absent in the ground state of the system, the impurity spin 
vanishes. Thus, as in the metallic Anderson system, the Kondo effect takes place at an arbitrary particle-impurity 
coupling. 

Solving Eq. (11), we immediately find n d — 1. Thus, at e d < —T 2 /4f3 the impurity level is entirely filled out and its 
occupancy does not depended on a position of the impurity energy with respect to the Fermi level. Thus, the behavior 
of the impurity occupancy in the gapless host is not described by a universal function of the renormalized impurity 
energy e* d but it depends essentially on the bare parameters of the model e d and f3. Making use of the terminology 
of the Anderson model, we will call this regime with the entirely filled impurity level the Kondo (or local-moment) 
regime, despite it disappears in the limit (3 — > 0, where our model reduces to the conventional Anderson model. 
However, for quite large 5, and even at 5 < 1, the Kondo regime describes the system's behavior almost at all negative 
td-, except a very narrow region near the Fermi level, where < b < j. 

The region < b < j corresponds in our case to the mixed-valence regime, where the impurity occupancy is 
changed from n d ss in the empty-level regime at e d > to n d = 1 precisely in the Kondo regime at e< — T 2 /4/3. The 
bare energy of charge complexes £o(A) = 4ra;(A) in this regime is negative at all A, except the point A = —5, where 
x(A) = 0. The renormalized energy of charge complexes at zero temperature is found from the thermodynamic BA 
equation 

/•Qi roc 

£(A) = 4r^(A) - / dX'a(X — A')£(A') - / dX'a{X - A')C(A'), (13) 



J — 00 J Qi 

where Qi and Q2 are defined as the zeroes of £(A), £(Qi) = £(Q2) = 0. This equation, as well as a corresponding 
equation for the density of states er(A), is hardly solved analytically, and a numerical analysis is required. 

In the metallic limit, f3 — > 0, the lower edge of the mixed- valence regime in the gapless system is shifted to 
—00. Correspondingly, the mixed- valence regime of the gapless system is extended to the conventional mixed- valence 
(0 < rid < 1) an d Kondo (rid ~ 1) regimes of the Anderson impurity in a metallic host. 

In summary, using the BA we have studied the ground state properties of an U — > 00 Anderson impurity embedded 
in a gapless host with an energy independent hybridization and the density of band states given in Eq. (1) with the 
power r — 1. As in the metallic version, the ground state of the system has been shown to be composed of charge 
complexes only. Since each complex contains two charge excitation and one spin wave, the total spin of a complex 
equals zero. Therefore, at zero temperature the impurity spin vanishes, and the Kondo effect takes place at arbitrary 
parameters of the model. 

However, the appearance of extra energy parameter /3, characterizing the size of region with an unconventional 
behavior of the density of band electron states leads to significant reconstruction of the density of states of charge 
complexes in the ground state of the system. This reconstruction results in significant changes in the behavior of the 
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impurity occupancy as a function of the bare impurity level energy. The empty- level (rid ~ 0) and Kondo (rid = 1) 
regimes are extended to almost all positive and negative values of ea, respectively. While the mixed- valence regime 
(0 < rid < 1) is squeezed to a narrow region — r 2 /4/3 < < 0. 

Only in the empty-level regime the impurity occupancy is a universal function of the renormalized impurity energy 
(*d = e d + £ l n Y + A which contains the parameter /?. The behavior of rid in the mixed- valence and Kondo regimes 
is not universal that demonstrates nonuniversal properties of the system. 

The latter explains why the Kondo screening is absent - or exists only at quite a large electron-impurity coupling 
- in earlier studies Jl|-01 based on scaling arguments. Nevertheless, by using both poor-man's scaling and numerical 
RG calculations, Gonzalez-Buxton and Ingersent have derived the behavior of the impurity occupancy which is 
qualitatively close to the picture described above. However, in their studies the particle-hole symmetry is a key 
factor in determining the low-temperature physics of the system, while in the BA solution, the ground state properties 
are obviously insensitive to this symmetry. 

In the BA analysis, the power r must be fixed. The spectrum of the system is determined by the functions k(u>) and 
h(u>), and it is essentially different at different r. At r ^ 1, the spectrum of Bethe excitations is enriched that could 
result in qualitatively novel physical properties of the system. Although it seems very difficult, if not impossible, to 
propose any scenario of destroying the Kondo screening of an integrable Anderson impurity. 

I thank S. John for stimulating discussions. 
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